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ABSTRACT 

The nearby (distance~350-400 pc), rich Vela OB2 association, includes y 2 Velorum, one of the most massive binaries in the solar 
neighbourhood and an excellent laboratory for investigating the formation and early evolution of young clusters. Recent Gaia-ESO 
survey observations have led to the discovery of two kinematically distinct populations in the young (10-15 Myr) cluster immediately 
surrounding y 2 Velorum. Here we analyse the results of Gaia-ESO survey observations of NGC 2547, a 35 Myr cluster located two 
degrees south of y 2 Velorum. The radial velocity distribution of lithium-rich pre-main sequence stars shows a secondary population 
that is kinematically distinct from and younger than NGC 2547. The radial velocities, lithium absorption lines, and the positions in 
a colour-magnitude diagram of this secondary population are consistent with those of one of the components discovered around y 2 
Velorum. This result shows that there is a young, low-mass stellar population spread over at least several square degrees in the Vela 
OB2 association. This population could have originally been part of a cluster around y 2 Velorum that expanded after gas expulsion or 
formed in a less dense environment that is spread over the whole Vela OB2 region. 

Key words. Stars: formation - Stars: pre-main sequence - Techniques: spectroscopic - open clusters and association: Vela OB2 - 
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1. Introduction 


Observations of star forming regions (SFRs) in the solar neigh¬ 
bourhood show that very young stars can be found in both high- 
density, gravitationally-b ound clusters and low-density, loos e 
stellar associations fe.g.. ICarpenteil 120001 : lLada & Lada| 2003). 
The origin of these young stellar populations is still being 
debated. All stars could form in massive dense clusters that 
expand to a larger volume after gas expulsion due to stel- 
lar winds, radiation pressure, and supernova explosio ns (e.g., 
Goodwin & Bastianl 120061 : iBaumgardt & Kroupal 120071) . Alter¬ 
natively, stars may form in a wide range of environments, 
from the dense core s of massive cluster s to sparse associations 
llBressert et alJl201fl : Wright et alJl2014h . High-resolution spec¬ 
troscopic surveys of SFRs and young clusters are powerful tools 
for investigating this topic, since stellar radial velocities (RVs) 
allow us to characterize their current dynamical state, establish- 
ing whether a SFR is bound or expan ding into the field (e.g., 
ICottaar et al.ll2012al : lJeffries et al.ll2014l) . and for identifying and 
disen tangling multiple populations seen in the s ame line of sight 
(e.g. JJeffries et al.ll20Q6bfc A1 ves & Bou\l20l2l) . 

The Vela OB2 association consists of 93 early-typejitars 
selected using Hipparcos proper motions by Ide Zeeuw et alj 
( 1999 ). which are distributed over an area of ~ 180 deg 2 . The 


* Based on observations made with the ESO/VLT, at Paranal Ob¬ 
servatory, under program 188.B-3002 (The Gaia-ESO Public Spectro¬ 
scopic Survey) 

** Table 1 is only available at CDS via anony¬ 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 


mean Hipparcos distance to the Vela OB2 stars is 410+12 
pc, but its most massive star, y 2 Velorum, is closer accord¬ 
ing to two interferometric observations of its orbit (336^ pc, 
iNorth et all[20071 368^f pc, iMillour et al.ll2007l) . The associa¬ 
tion is surrounded by a dust shell of radius ~7.5 deg, centred 
on y 2 Velorum, which is probably powering its expansion ( Sahui 
119921: lPetterssonll2008h . Recently, the young cluster (hereafter, 
the Gamma Vel clust er) around y 2 Vel has been observed by 
the Gaia-ESO Survey (iGilmore et al.ll2012l : iRandich & Gilmora 
120131) . The results of these observations, covering about one 
square degre e around the massive binary, are discussed in 
[Jeffries et al.l (120141) . They found that the cluster is composed 
of two kinematically distinct populations, which are older than 
previously thought (>10 Myr) and have an age difference of ~1- 
2 Myr. The older population (Gamma Vel A) is more concen¬ 
trated around y 2 Vel and is roughly in virial equilibrium, while 
the younger population (Gamma Vel B) is extended, clearly un¬ 
bound, and supervirial; namely, the total kinetic energy of the 
cluster is higher than in a bound system in equilibrium. 


The cluster NGC 2547 (RA = 08hl0m, DEC=-49dl2m) is 
located about two degrees south of y 2 Vel (~ 10 pc at the dis¬ 
tance of y 2 Vel). Several distan ce estimates, from 360-480 pc , 
are available in the literature (Ida Ll 1 19821 IShohhrookl 1 1 9861 
iLvra et all 1 20061 iNavlor & Jeffriesl 20061 van LeeuwetJ 20091) . 
An accurate age of 3 5+3 Myr has been derived from the lithium- 
depletion boundary (iJeffries & Oliveira! 12005 ). The Gaia-ESO 
cluster targets are selected to well sample the age-distance- 
position parameter space. NGC2547 is particularly important 
given its age (one of the few clusters with an age between 10 
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and 50 Myr) and relatively short distance to the Sun. Here, we 
report the discovery of a secondary population that is younger 
and kinematically distinct from NGC 2547. We show that this 
second component has a similar age to Gamma Vel B and is 
kinematically indistinguishable from it, providing evidence for a 
very extended low-mass component of the Vela OB2 association. 



Fig. 1. RV distribution for all the stars with EW(Li)>100 mA. The three 
lines describe the best-fit models obtained using: (a) one stellar popula¬ 
tion with a Gaussian distribution (green dotted line), (b) one population 
with a Gaussian distribution and a second one with the same distribu¬ 
tion of Gamma Velorum B (blue dashed line), and (c) two populations 
with Gaussian distributions (red continuous line). Stars with > 80 per 
cent probability of belonging to population B according to model (b) are 
highlighted with blue shading. The vertical dotted black lines indicate 
the centr al velocity of Gamm a Vel A (left) and Gamma Vel B (right) 
found bv lJeffries et alJ d2014f) . 


2. Observations and data 

The Gaia-ESO survey ob servations (see iGilmore et alJ 120121 : 
iRandich & Gilmore! 12013i) were carried out with the multi¬ 
object optical spectrograph FLAMES at the VLT, composed of 
a robotic fibre positioner feeding the GIRAFFE (R~ 17,000) and 
UVES (R~47,000) sp ectrographs, with 132 and 8 fibres respec¬ 
tively dPasquini et alj2003) . 

In this paper, we only use GIRAFFE data. The GIRAFFE 
candidate members are selected using a common strategy dis¬ 
cussed in Bragaglia et al. (in prep.). Specifically, for NGC 2547 
we used a two-step process: 1) we identified the regions includ¬ 
ing all the known cluster members in optical colour-magnitude 
diagr ams (CMDs), using optical photometr ic data in the litera¬ 
ture (iNavlor et al.ll2002tljeffries et alJl2004l) and a catalogue of 
members based on X-ray observations dJeffries et alj|2006al) : 2) 
we selected all stars in these regions with 11 < V < 19 (or 
10.5 < I < 16 in the absence of V photometry). The selected 
catalogue of candidate cluster members includes 467 stars dis¬ 
tributed over an area of about 1 deg 2 . We observed 450 candi¬ 
dates in 16 FLAMES fields (diameter ~25 arcmin each) with GI¬ 
RAFFE, during two separate runs in January and February 2013. 
Many stars were included in two or more overlapping fields and 
observed multiple times. In each field at least 20 fibres were used 
to acquire sky background spectra. All the observations were 
performed with the HR15N filter covering 647-679 nm. Expo¬ 
sure times were 1200 or 3000 s, depending on the brightness of 
the targets. 

The survey spectra were processed and analysed by 20 
working groups organized in a workflow, that is described in 
iGilmore et al.l d2012l) . GIRAFFE spectra were reduced with a 
pipeline developed at the Cambridge Astronomical Survey Unit, 


which also calculates radial velocities (RVs) and projected ro¬ 
tation velocities. A more detailed description of the data reduc¬ 
tion and an asses_sment of the accuracy and the precision of the 
RVs is given in iJeffries et al . (2014 ) and Lewis et al. (in prep.). 
The reduced spectra of young stars were analysed by a dedicated 
working group composed of several teams, which derive the stel¬ 
lar parameters (effective temperatures, gravities, chemical abun¬ 
dances, accretion tracers, as well as equivalent widths of the 
Li i 670.8 nm line - EW(Li)). The final parameters are weighted 
means of the results provided by the teams after outliers are re¬ 
jected. A detailed description of the methods used and_a dis¬ 
cussion of accuracy and precision are given by Lanzafam e et al.l 
( 2015 ). The data presented in this paper are part of the second 
internal data release (GESviDR2) of July 2014 and are reported 
in Table 1, which is available at the CDS. 
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Fig. 2. Equivalent width of Li 670.8 nm absorption as a function of R-I 
(corrected for the extinction E(B-V) = 0.06 derived bv IClari 1E9H) 
for stars with > 80 per cent probability of belonging to population A 
(empty diamonds) or population B (red circles). The black dashed lines 
are isochrones of Li depletion calculated from the lBaraffe et al.l ( fl998i) 
models at 10, 32, and 200 Myr, and the horizontal dotted blue line in¬ 
dicates the threshold used to select the pre-main sequence stars. The 
x-axis on the top of the panel rep orts a sc ale of effective temperature de¬ 
rived from R-I using the lPecaut & Mamaieki d201 3l ) colour-temperature 
transformations. 


3. Population B of NGC 2547 

The strategy for the selection of targets for young clusters in the 
Gaia-ESO Survey is tailored to observe an unbiased and nearly 
complete sample, but can lead to significant contamination by 
field stars (e.g.. IJeffries et al.ll2014l) . To cleanly study the kine¬ 
matics of the young NGC 2547 population and search for ad¬ 
ditional young populations, we need to remove these contami¬ 
nants. 

Since pre-main sequence (PMS) stars with T e ff ~4000 K 
completely deplete their photospher ic lithium in <100 Myr 
(ISiess et al.ll2000t iBaraffe et al.lll998h . a filtered sample can be 
created by excluding all stars with no evidence of Li absorption 
at 670.8 nm. This approach isolates a sample of young stars at 
cool temperatures, but will not select all members of NGC 2547, 
since many M-type PM S stars at an age of 35 M yr will also have 
depleted their lithium ( IJeffries & Qliveirall2005l) . 

Figure [I] shows the RV distribution of all the observed stars 
with EW(Li)>100 mA. The distribution is characterized by two 
peaks at ~13 km s _I and ~20 km s -1 , respectively. The position 
of the former and most prominent peak is consistent with the 
meanRV ofNGC2547 (12.8+1 kms ') derived bv lJeffries et al.l 
( 2000 ). while none of the previous spectroscopic studies of this 
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Table 1. Results from maximum likelihood modelling of the radial velocities. 



One component 

One component+y Vel B 

Two components 

RV a (km s 1 ) 

14.68 + 0.44(14.79) 

12.94 + 0.17(12.87) 

13.05 + 0.21 (13.13) 

cr A (km s~') 

3.07 + 0.36 (2.95) 

0.83 +0.17(0.74) 

0.97 ± 0.23 (0.95) 

RV b (km s' 1 ) 


18.88 

18.9 + 0.95 (19.68) 

<tb (km s~') 


1.6 

1.83 ± 1.40 (0.61) 

h 

1.00 

0.69 ± 0.06 (0.70) 

0.72 ± 0.09 (0.75) 

In (Lmax) 

-194.9 

-174.6 

-173.81 


Notes. Symmetric 68 per cent confidence intervals for one parameter of interest, the value at the maximum likelihood is given in parentheses. 
The parameters RV A , RV b , cr A and cr B represent the mean and the standard deviation of the normal distributions describing the intrinsic RV 
distributions of the single stars in population A and B. The parameter f A is the fraction of stars belonging to population A, while In L max is the 
maximum log-likelihood of the fits. 
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Fig. 3. Comparison between the multiple populations discovered in 
NGC 2547 and Gamma Vel B. The top panel shows the EW(Li) 
as a function of V-I, while the bottom panel shows a V vs 
V-I CMD. Photometry has been corrected according to the extinc¬ 
tion towards NGC2547 (E(B-V)=0.06) and Gamma Velorum (E(B- 
V)=0.055). The x-axis on the top reports an effective temperature scale 
derived from V-I, using the colour-temperature transformations de¬ 
scribed in iPecaut & Mamaiekl J2013 I). In both panels empty squares, 
filled red circles, and blue circles indicate the stars of Gamma Vel B, 
NGC 2547 B, and NGC 2547 A respectively. EW(Li) upper limits in 
NGC 2547 A are reported with downward triangles. Only a subsample 
of stars of NGC 2547 is shown, since V band photometry is not avail¬ 
able for all targets (see footnote on page 4). 


cluster mention a kinematically distinct secondary population 
with very_similar velocity to Gamma Vel B (18.9 ± 0.5 km s 1 
from lJeffries et alJi2014l) . To test the hypothesis that the popula¬ 
tion associated with the second RV peak (called NGC2547 B in 
the rest of the paper) and Gamma Vel B have the same origin, we 
model the RV distribution shown in Fig.Q]using the same maxi¬ 
mum likelihood technique as described in detail in ljeffries et ali 


(2014j), which is based on a method developed by ICottaar et ali 
(l2012bi) . We assume that the intrinsic RV distribution of a sin¬ 
gle population is described by a Gaussian that is broadened by 
binary motions and the uncertainties in RV measurements. The 
RV uncertainties are calculated using equation 1 in ljeffries et al.l 
(12014 . since NGC 2547 and Gamma Vel were observed with the 
same instrument and processed with the same pipeline. We tried 
three different models: a) a single population with mean veloc¬ 
ity RVa and dispersion <r A ; b) two populations, one with mean 
velocity RVa, dispersion cr A , and fraction of the total population 
/a as free parameters, and a second one with a fixed mean ve¬ 
locity (RVb = 18.88 km s -1 ) and dispersion (<xg = 1.6 k m s 1 ), 
which are the values found for Gamma Vel B by iJeffries et al.l 
(12011 : c) two populations with five free parameters (mean ve¬ 
locities RVa and RVb, dispersions cr A and cr B , and ratio between 
the two populations f A ). The most likely distributions derived 
by the three fits are shown in Fig. Q] The best values of the 
free parameters and the maximum log-likelihood are reported 
in Table Q] Since the models are nested, we can perform a like¬ 
lihood ratio test to establish which model describes the RV dis¬ 
tribution better. The test rejects model (a) when compared to the 
others with a probability > 99 per cent, while the difference in 
maximum log-likelihood of models (b) and (c) is not significant 
(P(ln L maxA - In L maX B) > 46 per cent). We conclude that the RV 
distribution of the Li-rich members of our sample is composed 
of two kinematically distinct populations and that the RV dis¬ 
tributions of NGC 2547 B and Gamma Vel B are consistent, in 
agreement with the initial hypothesis. 

The cluster NGC 2547 (age~35 Myr) is significantly older 
than Gamma Vel (age-10-15 Myr), therefore we expect to see 
significant differences between the EW(Li) measured for pop¬ 
ulations A and B of NGC 2547 if the latter is associated with 
Gamma Vel. In Fig. [2] we plot the EW(Li) as a function of R-1, 
together with isochrones at 10, 32, and 200 Myr derived from 
the lBaraffe et al.l (119981) models. Different symbols are used for 
stars that - according to their RV - are associated with one of the 
two populations with > 80 percent confidence. As expected, in 
the colour range 1.0 < R — I < 1.8, we observed very few stars 
belonging to NGC 2547 A, because at an age of 35 Myr they are 
Li-depleted, and so were excluded from our initial sample selec¬ 
tion. However, we observed several stars belonging to NGC 2547 
B in this colour range, proving that this population is younger 
than 35 Myr. The evidence of this age difference is supported 
by the presence in NGC 2547 B of a star (2MASSJ08104437- 
4939001) classified as a strong accretor on the basis of its Her 
emission (EW(Har)=78 A and a line width at 10 per cent of 440 
km s _1 ). Such strong accretors are extremely rare in a 35 Myr 
old cluster, and there are none in NGC 2547 A. 
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In Fig. [3 we compare the two NGC 2547 populations (stars 
of NGC 2547 A have been selected on the basis of their RV, in¬ 
cluding Li-depleted stars) with Gamma Vel m The top panel 
shows that the EW(Li) measured in Gamma Vel B and NGC 
2547 B are consistent, suggesting that they are of similar age and 
younger than NGC 2547 A, because Li depletion occurs rapidly 
at V - I ~ 2.5. The bottom panel shows that Gamma Vel B and 
NGC 2547 B share the same locus in the CMD, which, assum¬ 
ing that they are coeval, suggests they are at a similar distance. 
The exceptions are two stars in the upper left-hand comer of 
the CMD, that match the sequence of NGC 2547 A. However, 
Li absorption cannot reliably select young (< 100 Myr) stars at 
T eS >5000 K (V - I <1), so these two objects, which have RV 
(21.65 + 0.19 and 23.10 ± 0.25 km s~') might be contaminating 
field stars or are possibly binary members of NGC 2547. 

4. Discussion and conclusions 

Ijeffries et al.l ( 2014 1 propose several scenarios to explain the 
presence of two populations around y 2 Velorum, concluding that 
population A is the remnant of a bound cluster formed around the 
massive binary, and population B is a dispersed population from 
the wider Vela OB2 association. We have discovered 15 stars that 
appear to belong to Gamma Vel B and are located 2 deg (~10 pc, 
assuming a distance of ~340 pc) south of y 2 Velorum, demon¬ 
strating that_Gamma Vel B extends far beyond the area studied 
bv Ijeffries et al.l (2 0141 ). Two scenarios might explain this result: 
a) the young low-mass population observed towards NGC 2547 
was part of a denser cluster around y 2 Velorum, which - after 
the formation of the massive binary and the expulsion of residual 
gas - expanded into a larger volume; b) these stars were born in 
a low-density diffuse environment, which also formed the other 
members of the Vela OB2 association. 

According to the relation between cluster mas s and the mass 
of its most massive star proposed by [Weidner et al. feoidl . 
y 2 Velorum should be surrounded by a cluster with a mass 
(~ 1000 Mg), which is much higher than the present total mass 
(~ 100 Mg) estimated bv lleffries et al.l (i2014l) . It therefore seems 
likely that most of the original cluster that formed around the 
massive binary has expanded into a larger volume. Assum¬ 
ing that the cluster started to expand soon after the formation 
of the massive binary (age~5.5 ± 1 Myr, lEldri dgei i2009i) and 
considering the intrinsic velocity dispersion of Gamma Vel B 
(~1.6 km s-'), it is possible that the fastest stars formed around 
y 2 Velorum have since moved ~10 pc. This possibility has been 
confirmed by recent N-body simulations aimed at investigating 
the origin of the two kine mati cally distinct populations discov¬ 
ered around y 2 Velorum dMapelli et al.ll2015l) . Specifically, these 
simulations show that Gamma Vel B is supervirial and extends 
up to a projected distance of ~10 pc from the massive binary. 
The predicted stellar density in the outer region of the cluster 
is consistent with the number of stars found in NGC 2547 B. 
Furthermore, NGC 2547 B is only composed of low-mass stars 
(R-I > 1.2, i.e., below ~ 0.5 M Q ). This supports the cluster 
expansion hypothesis, since in this scenario we should observe 
a higher concentration of stars in the cluster centre and mass 
segr egati on has been observed in many massive clusters (e.g., 
Ide Griis et al.ll2002lf . 

However, not all the evidence fits this first scenario. The pres¬ 
ence of Li-depleted stars proves that Gamma Vel is older (age > 

1 A complete photometric catalogue in the V, R and / bands is not 
available for all targets in Gamma Vel or NGC 2547. Therefore not all 
the stars plotted in Fig.[2]are included in Fig.[3]and vice-versa. 


10 Myr) than the massive binary. The putative expanding clus¬ 
ter defined by Gamma Vel B is also offset by 2 km s _1 in RV 
from Gamma Vel A, which appears more centrally concentrated 
around y 2 Velorum. Furthermore, the other early-type members 
of the Vela OB2 association are spread over an area of radius ~7- 
8 degrees (~50 pc at a distance of 400 pc), and considering the 
timescale of cluster expansion, it is unlikely that they all formed 
in a much smaller region and spread out after the formation of 
the massive binary. Therefore, the second scenario seems plau¬ 
sible, but fails to explain why only very low-mass members of 
Gamma Vel B have been found towards NGC 2547 - the contrast 
in the distribution of Gamma Vel B and NGC 2547 B members 
in the bottom panel of Fig. 2 is striking. 

The Gaia satellite will soon provide accurate parallaxes and 
proper motions dde Bruiinel2012l) for the whole Vela OB2 region 
to V ~ 19. Therefore, it will be possible to perform an unbiased 
census of the young stellar populations and study their three- 
dimensional spatial and kinematic structure. This will allow us 
to answer this and other questions concerning the star formation 
process in Vela OB2 uncovered by the Gaia-ESO Survey obser¬ 
vations. 
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